Excitation modes of vortices in sub-micron magnetic disks 



R. Zarzuelai'2, E. M. Chudnovsky^, J. Tejada^ 
^ Departament de Fisica Fonamental, Facultat de Fisica, 
Umversitat de Barcelona, Avinguda Diagonal 645, 08028 Barcelona, Spam 
^Physics Department, Lehman College, The City University of New York, 
250 Bedford Park Boulevard West, Bronx, NY 10468-1589, U.S.A. 
(Dated: September 7, 2012) 

Classical and quantum theory of spin waves in the vortex state of a mesoscopic sub-micron mag- 
netic disk has been developed with account of the finite mass density of the vortex. Oscillations of 
the vortex core resemble oscillations of a charged string in a potential well in the presence of the 
magnetic field. Conventional gyroscopic frequency appears as a gap in the spectrum of spin waves 
of the vortex. The mass of the vortex has been computed that agrees with experimental findings. 
Finite vortex mass generates a high-frequency branch of spin waves. Effects of the external magnetic 
field and dissipation have been addressed. 
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I. INTRODUCTION 

Recent advances in optical and electron-beam lithog- 
raphy offered possibility to fabricate arrays of micron 
and submicron-size magnetic structures with controlled 
magnetic properties. Among such structures are meso- 
scopic circular disks of soft ferromagnetic materials. Ar- 
rays of such disks, as well as individual disks, have 
been intensively studied -'^^^ due to their unusual magnetic 

properties and potential for technological and biomedical 
application^sEol, 

Micron-size circular disks exhibit a wide variety of 
magnetic equilibrium configurations due to geometric 
constraints on the spin fielcP^. Their applications are 
based on static and dynamic properties of one of the es- 
sentially non-uniform ground states, the vortex state. It 
is characterized by the curling of the magnetization in the 
plane of the disk, leaving virtually no magnetic "charges". 
The very weak uncompensated magnetic moment of the 
disk sticks out of a small area confined to the vortex core 
(VC). The diameter o f the core is comparable to the ma- 
terial exchange lengtlPE^. The low frequency dynamics 
of the vortex state is due to the gyrotropic mode, con- 
sisting of the spiral-like precessional motion of the VC as 
a whol^^^'^, and it is intrinsically distinct from conven- 
tional spin wave excitations. 

Because of the strong exchange interaction among the 
out-of-plane spins in the VC, it behaves as an indepen- 
dent entity. The research on excitation modes of vortices 
has focused on the low-frequency gyroscopic mode that 
describes circular motion of the vortex about the center 
of the disk. It can also be viewed as the uniform pre- 
cession of the magnetic moment of the disk due to the 
vortex. The natural question is whether the gyroscopic 
mode allows spatial dispersion similar to spin waves of 
finite wavelength in ferromagnets. The aim of this paper 
is to study spin waves related to the gyroscopic motion 
of the vortex. Such a wave is shown in Fig. 1. It must 
exist due to finite elasticity of the vortex provided by the 
exchange interaction. 




Figure 1: Gyroscopic spin wave in the vortex state of a meso- 
scopic magnetic disk. 



Most of the research on the gyroscopic motion of vor- 
tices in circularly polarized disks ignores the inertial mass 
of the VC. Such a mass has a dynamical origin stemming 
from the variation of the shape of the VC as it moves in- 
side the disk. Meantime, experimental studies of vortex 
oscillations in micrometer permalloy rings hinted to- 
wards a non-negligible vortex mass of order lO^^^kg. On 
a theoretical side the vortex mass has been previously 
computed in a two-dimensional Heisenberg model with 
anisotropic exchange interactiorP^HH] jjj disks made of 
soft magnetic materials that have been experimented 
with, the exchange interaction is isotropic. We will show 
that in this case the finite mass density of the vortex 
originates from the geometrical confinement of the spin 
field and the magnetic dipole-dipole interactions. 

In this paper we will derive the generalized Thiele 
equation that describes spin waves in the vortex core of 
finite mass density and will obtain the spectrum of such 
waves. We will show that the conventional gyroscopic 
mode, ijQ, appears as a gap in the spectrum of the spin 
waves in the vortex, uj(q) — luq + aq^, when the vortex 
mass is neglected. From the mathematical point of view 
the above problem resembles the problem of the motion 
of a charged string in a potential well in the presence of 
the magnetic field. The latter problem is a generaliza- 
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tion of the problem of Landau levels of an electron in 
a two-dimensional potential well in the magnetic field. 
We will show that this problem has a nice exact solution 
for quantized oscillations of the string, thus providing 
the spectrum of magnons in the vortex in the quantum 
regime as well. Classical and quantum solutions for the 
spectrum of excitations of the vortex lead to the same 
dispersion law, Lo{q), in the limit of small q. 

The paper is structured as follows. In Sec. II a La- 
grangian formulation of the problem is presented. For- 
mal derivation of the massive elastic Thiele's equation 
that allows deformations of the vortex line is given in 
Section III. The spectrum of spin waves in the vortex 
core is obtained in Section IV. We show that a finite 
mass of the vortex results in the additional excitation 
mode that is absent in the case of zero mass. Quantum 
mechanical treatment of magnons in the vortex core is 
developed in Section V. The vortex mass in a circularly 
polarized disk is computed in Sec. VI and is shown to be 
in good agreement with experimental findings. The field 
dependence of the vortex excitation modes and effects of 
dissipation are discussed in Sec. VII. Sec. VIII contains 
final conclusions and suggestions for experiment. 



II. LAGRANGIAN MECHANICS OF THE 
VORTEX CORE 

We shall describe the vortex line by the vector field 
X = {x,y), where x{t,z) and y(t,z) are coordinates of 
the center of the vortex core in the XY plane. Landau- 
Lifshitz dynamics of the fixed-length magnetization vec- 
tor M(9,<i>) = Ms(cos$sin9,sin<I>sin9,cos8) follows 
from the LagrangiarP^l 



dz d^r 
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(A$) cose- 5(6, $,a,e,a,$) 



(1) 



where £{Q,<^,dzQ,dz^) is the energy density. The de- 
pendence of the Lagrangian on the partial derivatives 
of the angular coordinates comes from the elastic nature 
of the vortex core. It is contained in the total energy, 
£, that takes into account interaction between different 
layers of the vortex line, see below. 

The spatial dependence of angular coordinates (8, <i>) 
for the vortex state is given by = 8(t;r, z) = 8(r — 
X{t,z),t) and $ = $(t;f,z) = <^{f ~ X{t, z),t). We 
only consider long-wave solutions that do not deform the 
vortex core in any z-cross-section of the disk. This means 
that the angular coordinates depend on t and z via the 
coordinates of the vortex core X{t,z). The covariant 
derivative with respect to time, -Dt$, along the vortex 
core is given by 



X{t,z) 



-X{t,z) ■Vf.'S>{r~ X{t,z)) (2) 



where "dot" denotes partial derivative with respect to t. 

Taking all these considerations into account , the above 
Lagrangian becomes 



= / dzC 



C t]X,X,dzX 

with the Lagrangian density being 

M. 



(3) 



c 



= d 
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-X{t,z) ■ V.p$) cos 8 



f(8,$,9,8,a,$) 



(4) 



Thus the generalized momentum densities are given by 

6C 



t,z;X,X,dzX 



t,z;X,X,dzX 



5 (X(t,z)) 

/ dV(V,F$)cos8 (5) 

7 J 

SC 

d(^dzXit,z) 
6uj{X, dzX) 



S{dzX{t,z)) 



(6) 



with uj{X,dzX) 



dV £(8, $,9^8,92$) being the 



linear energy density. The dynamics of the vortex core is 
governed by the Euler-Lagrange equation, 



SC 



6X(t,z) 



(7) 



Notice that 



d 



DtA{t, z- r, v) = i{t, 2)(A) = ^A(f - X{t, z), X{t, z)) 

= -X{t, z) ■ VfA(f - X{t, z),X{t, z))+ 

X{t, z) ■ VffA(f - X{t, z),Xit, z)) (8) 

means covariant derivative along the curve that is tangent 

to the vortex core, ^(t, z). All terms involving X{t, z) and 

z) in the Euler-Lagrange equation come from Dtllt, 
which is given by 



An* 



7 
7 



dVA (cos8Vf.$) 



d^f Vf 



COS 8 



—Xjdj cos 8 + Xjdj cos 8 



Ci ^^ij Xj Ci Kij Xj , 



(9) 
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where 



and di = 



7 



dV 



{didj^) cos e + cos e 

(9^9^$) cos e + {di^)dj cos 9 



(10) 



We want CiKijXj to be of the form pc 'x X which re- 
sults in the identity SijkPG.j — Kik- From this we obtain 

Pgj — —-^ikjKik, which translates into the vector form 



PG 



27 
27 



dVe 



dV 



^kJ cos e + cos 6) 

(V,-^- X Vf?$) cos e+ 

V,?<I> X cos6 



(11) 



To compute the mass tensor and the gyrovector we 
have to find the solutions (0, $) of the Landau-Lifshitz 
equation in the low dynamics regime that is characterized 

by the condition |X| ^ 1. In this regime solutions can 
be expanded as a perturbative series on the differential 

speed, \X\, of the vortex core 

e{t, z;r)^ e(°) (z; r) + e^^) {t,z;r) + ... 

<i>{t,z;r) ==$(°)(z;r) + $(i)(i,z;r) + ... (12) 

Notice that the zero-th order is time independent, other- 
wise the gyrovector would depend on time. 

The approach that neglects deformation of the vortex 
core in any z-cross-section of the disk is correct only for 
weak deviations of the centerline of the vortex core from 
the straight line along the Z-axis. We now proceed to 
the study of the Landau-Lifshitz equation for the set of 
variables (0, $) in such weak bending regime. It can 
be obtained by applying the variational principle to the 
Lagrangian density 



M, 



(A$)cose 



£(e,$,v,-e,v,-$,a,e,<9,$) 



(13) 



Notice that 



£■(6,$, V^.e, Vf.$,9^e,9^$) = £'xy(e, $, V^^O, V^.$) 

+ £,iiQ,^,d,e,d,^) (14) 

with V,-5-0, Vp$) being the sum of the ex- 

change, anisotropy and dipolar energy responsible for the 
formation of the vortex, and 

(15) 



being the elastic energy in which A^/f is a constant. It 
describes contribution of the exchange and dipolar forces 
to the elasticity of the vortex line, with the exchange 
playing a dominant role. Consequently, with good accu- 
racy, Aeff can be identified with the exchange constant 
A. 

The set of dynamical equations for (8, $) is 



D 



SC 



6C 



2Aeff [sin2ea,e<9,$ + sin^ e9>] 
Ms dcos 9 S£xY 
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dt 



(5$ 



and 




5C^ 
56 







= (16) 



= 



-2A 



Performing a Fourier transform 
1 



$(t,z;rO 
e{t,z;r) 



/2n 
1 

/2^ 



dq $(i, q; f)e*«^ 
dqQ{t,q;'r)e"" (18) 



we obtain the following set of equations for the pair 
(6,$): 



A, 



sin 26 7k- 9 ★ <i> -t- sin 9 ★ sin 6 ★ $ 



Ms d cos 9 6£xY 



7 



dt 



5<P 



(19) 



1 



, -sin 9* 

/27r 7 dt 



d<P 5£ 



XY 



5Q 



+ 



2A, 



6 



sin 26 
47r 







(20) 



where ★ means Fourier convolution. 

For a small bending of the vortex core, the boundary 
conditions on the angle 6 are the same as in the rigid VC 
case, i.e.: 6 ~ or tt in the limit f <^ Aq and 6 ~ tt/2 in 
the limit f 3> Ap , with Aq = y^A/M^ being the exchange 

length of the material and where f = ||r — X(t, z)||2 is 
the radial distance from the VC center at any height z. 



Considering this two limits Eq. ( 19 ) becomes 



• Limit r <C Aq. In this case, sin 6 ~ and thus 
sin 269,992$ + sin^ 99^$ ~ 0. So we have the 
following equation in the Fourier space 



Ms d cos 9 
^ Jt 



S£ 



XY 

(5$ 



= 



(21) 
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Limit f ^ Aq. In this case, sinO ~ 1 and thus 
sm2&d^ed^<^ + sin^ 69^$ ~ ^2^. So we have the 
equation 



Ms d cos 8 



7 



dt 



5£ 



XY 



(5$ 







(22) 



Notice that in both Hniits sin 20 
Eq. ( 20 I becomes 



and so d^Q — 
Consequently, in the Fourier space 



1 AL 



/2n 



-sin 8 ★ 
7 dt 



5E 



XY 



(58 



(23) 



Finally, in the hmit of weak bending {Aef/q^ ^ 1), 
we can neglect the terms of the form 2Aeffq^£, in the 
above equations. In doing so, we recover the standard 
Landau-Lifshitz equations for (8, $) at any z layer, with 
the VC center depending on the value of z. Introducing 
now the perturbative series ( |18[ ) into the Landau-Lifshitz 

equation and splitting it into 0(|X|") terms, we obtain 
the equations of motion for the $(")/©(") terms. In the 
case of the zero-th and first order terms, we recover the 
static solution and the first perturbative solution for the 
rigid vortex (see Section VI). For the particular case of 
the zero-th order we obtain 

^o{x,y) = n^ta.ii'^{y - yy/x - Xy) 

c»6.(,-, '^*^» (24) 

where = ±1 is the vorticity of the magnetization of 
the disk and Ci , C2 are constants that can be obtained by 
imposing the smoothness condition on cos 80 at f = Aq 
up to its first derivative. The corresponding values are 
Ci = I and C2 = |pe. From all this we straightforwardly 
deduce that 

X V?$o = 2T:nyS'-^^ (^r - X(t, z)^ 

Vl^o = (25) 



III. ELASTIC THIELE EQUATION 

We now proceed to the computation of the gyrovector 
and the mass density tensor. Using (25 1 we obtain that 
the first term of Eq. (Ill equals {TmypMs/^)ez, where 



p = cos 8(0) = ±1 defines the direction of the polar- 
ization of the vortex core (8(()) = or tt). The second 
term is evaluated at the zero-th order of the perturba- 
tive expansion of the angular coordinates in the low dy- 



namics regime, Eq. (24 1, taking into account that in the 



weak bending regime the deformation of the vortex core is 
small and one can consider f ~ r because \ z) 1 12 <C 1. 
In doing so we obtain {TrnypMs/^)ez agairPS. Thus the 
gyrovector becomes pQ = pcpnyCz with 



PG = 27rMs/7 



(26) 



2 ln(i?/Ao), where R is the radius of the disk. Only 



Notice that pa is the gyrovector linear density as com- 
pared to the gyrovector in the Thiele equation for a rigid 
vortejP^I, 

Computation of the mass density tensor will be per- 
formed in Section VI. For circular polarized disks we 
show that this tensor reduces to a scalar, Mij = pM^ij, 
with the vortex core mass density given by pM — 

1 

47 

U!{X, dzX) contributes to the partial derivative SC/6X in 

the slow dynamics regime, because lit = PmX and so the 

term X ■ lit equals pmX^ ■ Consequently, the generalized 
Thiele equation becomes 

PMX{t,z)+X{t,z)-KpG + dznz+^x^^Q (27) 

The linear energy density uj{X , dzX) is the sum of the 
magnetostatic and exchange contributions in the z-cross- 
section, ujxYiX), and an elastic contribution due to the 
deformation of the vortex core line, ujei{dzX). Zeeman 
contribution will be considered later. The dependence 
on the vortex core co ordin ates on the ujxYiX) term for 
small displacements i^i^Ell 



ujxy{X) = ^pMOjlieoX'^, 



(28) 



where ujm = Pg/ Pm is the characteristic frequency of the 
system and eg = ljq/ujm is a dimensionless parameter. 
Recall that the conventional gyrofrequency ljq is defined 
asi^ 



ug 



J'xy{X^Q) 



PG 



(29) 



where /3 = L/R is the ratio of the thickness and the 
radius of the disk. This last expression is valid in the 
limit /? < 1. 

From the continuous spin-field model we know that 
uj,i{dzX) = A,ff J dV [(9,8)2 + sin2 8(a,$)2] (3^) 

Noticing that 9,8 = -VpQ-dzX and 9,$ = -Vp<^-dzX, 
and taking into account the vector identity, {A x B) ■ 
{C X D) ^ {A - C){B ■ D) - {A - D){B ■ C), with either 
A^ D = Vf8, B = C = -dzX OT A = D ^ Vf$, B ^ 
C — —dzX, we obtain the following relations 



(Vf8 • dzXy = {VfQf {^^^y ^ (^^-^ ^ 
(\7p'i> ■ dzXy = C^P^ f (dzXy - (^r^ X dzX 



(31) 



(32) 

The main contribution to the integral comes from the 
zero-th order in the perturbative expansion (18 1. Notice 
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that 60 (r) = Qo{r) and that V^^^o = ire,/,. As dis- 
cussed before, in the weak bending regime we can use 
approximation f ~ r. so that Vf<i>o x dzX — —^dzXriz 
and Vf9o x dzX = ^^dzX^iz, where Xr — Sr ■ X = 
X cos 9 + y sin 6 and X^ — e^ - X = —x sin 9 -\-y cos 9. The 
elastic energy density finally becomes 



UJel(dzX) ^A,ff / d' 



ax 

dz 



Aeff 

A, 



dr J \ dz 



2 ^ sin^ 60 f dXr 



dV 
dV 



\ dz 



=ttA, 



eff 



r dr 



\ 2 9 
d6o \ sin 9o 



dr 



dx_ 

dz 



(33) 



Finally, for the total energy density we obtain 

dzX) = ipM^l/eo^' + (39) 

and thus the generalized Thiele equation for an elastic 
vortex core line becomes 

PMX{t,z)-\dlX{t,z)+X{t,z)xpG+PM<-JMeaX{t,z) = 

(40) 



IV. SPIN WAVES IN THE VORTEX CORE 

Introducing the complex variable x — x — iy we can re- 
cast Eq. ( 40 1 as the following complex partial differential 
equation: 



PMX - A5^X + ipGX + PM^lieoX = (41) 



where the angular dependence of (dzXr)'^ and {dzX^)'^ 
has been integrated on 9. We can recast this energy 

density as ujei{dzX) = -A ; where A is the elastic 

constant given by 



A = 2t: Aeff I rdr 



dOo 
dr 



sin^ 60 



(34) 



Making use of the variable mo{r) — cos0o(r), we can 
rewrite the above equation as 



A = 2TTAeff 



• dr 



1 



/ drriQ 



1 — rriQ \ dr 



1 — TOq 



(35) 



and using the spatial dependence (24 1, we get 
2 



1 / dniQ 
1 — TOq \ dr 

2Cir2 



1 — TO 







2CirZJ_ 



4Ci 



r < Ao 



i^m^-l- = ^exp(-2r/Ao) + ^ r » Aq 



(36) 



Recalling that Ci = | and C2 = |pe and computing the 



integral (34 1 by splitting into two regions, [0, Ao] and 
[Aq, R], we obtain 



50 

A = 2^Ae//(--fln(i?/Ao; 



(37) 



Let Xo(^) be the equilibrium complex center of the 
straight vortex core line. In the presence of the wave 
it gets perturbed and becomes x{t, z) = Xo{z) +Xw{t, z), 
with llxiollz ^ llxolU- Switching to the Fourier trans- 
form, 



Xw{t,z) 



1 

2tt 



(42) 



we obtain the following equation for Xwii-^, <l)- 

[-PMUJ^ + - Pg^ + Pm^\i^q\ Xw[^; g) = (43) 

For non-zero amplitude of the wave the expression in 
the square parenthesis must vanish. This determines the 
spectrum of the waves: 



Pmw^ - Ag^ -t- pg^ - PMOjlieo = 



At pm 7^ one can normalize Eq. ( 44 1 to get 

A 



2" 9" 

PMl^M 







(44) 



(45) 



Solving this equation we obtain the spectrum of vortex 
core excitations: 



w±(g) 



'(l + 4eo) + 



4Ag2 
Pm^m 



± 1 



(46) 



In the weak bending regime we have A\q^ / pM'^\[ ^ 1 
and so we can expand the square root and obtain the 
following expression for the frequencies 



In the limit R ^ Aq the logarithmic term dominates and 
A becomes 



1 



A 



2 ^ ' VI -t- 4eo PG 



(47) 



A = 27r^e//ln(i?/Ao) 



(38) where we have used the relation pm^^m — Pg- 
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As will be shown in Section IV the parameter cq 
luq/uj]\i is normally small due to the smallness of /3 
L/R. Consequently 



« uJG + —q^ (48) 
Pg 

uj+iq) ~ UJM H (49) 

Pg 

With the help of Eqs. (|26| and (|38| with Ae// « A = 
M^A§ the above equations can be written in a transpar- 
ent form: 



a;_(q) 

^+iq) 



ujG+lMs{q^o?HR/^o) (50) 
i^M+lAUqAoflniR/Ao) (51) 



Note that the weak bending regime corresponds to 
qAo < 1. 



V. QUANTUM MECHANICS OF THE 
EXCITATIONS IN THE VORTEX CORE 

In this section we will show that excitations of the 
vortex core can be also obtained in a rather non-trivial 
way from the quantum theory as well. This problem is 
interesting on its own as it turns out to be equivalent to 
the problem of quantum excitations of a charged string 
confined in a parabolic potential and subjected to the 
magnetic field. 

It is straightforward to prove that the generalized 
Thiele equation (40 1 is the Euler-Lagrange equation as- 



sociated with the following effective Lagrangian density 
that can be derived from Eq. Q 

C{t, z; X, X, d,X) = ipM^' +X-Ap^- w(X, d,X) 

(52) 

where Ap^ is the gyrovector potential satisfying x 
Ap^ — -^pg- Thus the total Lagrangian becomes 



AzC 



Az 



1 



-PmX^ +X-Ap 



L0{X,d,X 

(53) 



Noticing that Wn{z)}neT<i 



271- 




sin (g„z) 



with 



g„ = —n, is a Hilbert basis of the function subspace 
ij 

W = {ip ^ £2(0, L), (p(0) = ip{L) = 0}, we can expand 
X as 



X{t,z) = Xoit) +J2Xnit)iPn{z), 



(54) 



where A'o(t) is the center of the undisturbed vortex and 

Xnit) =< X{t, z), (pn{z) >£2(o,L)- Introducing this ex- 
pansion in Eq. ( 53 I and taking into account the orthonor- 



obtain the following identity: 

citAxA ,fi„} \ 

\ L J nez+ L J riez+ / 



2 

E 

n>0 



-MXi 



Xo- Ao-\MujlieoX^ 



■^PM^n + Xn ■ An - -pA/CJ^.^eo^^ - -XqlX^ 



(55) 



where Z+ — {0} U N. M = PmL is the total mass of 
the rigid vortex line and An is the gyrovector poten- 
tial associated to the n-th coordinate X„, which satisfies 
V X Aq = —G and x A„ = — pb, n > 0, with 
G = pgL being the gyrovector of the rigid vortex. 

Applying the Laguerre transformation to the above 
Lagrangian we obtain the following expression for the 
Hamiltonian 



±-^{n^ - A^f + ]^Mul,eoXl 



E 

ri>0 



1 



2p 



M 



(n. 



An)- + ^Pm^^m eo + 



2 \ x?2 

(56) 



where Hq = MJ?o + a-nd n„ = pMXn+ An, n > are 
the corresponding canonical momenta. Notice that Eq. 
(56 1 shows that H splits into the direct sum ffimgz+^mj 



with Tim being the Hamiltonian defined over the phase 
space (XmjIIm). It has a structure of the form 

n' = ^^{Il-Af + lvu^Ux', (57) 

where {X, 11) are the canonically conjugate variables, 77, ^ 
are constants and A is the gyrovector satisfying V x A = 
— Xz, with X being a constant. It is important to point 
out that x/77 — lom in all cases. 

From now on we consider the case of the vortex core 
of a non-zero mass, (77 0) . It is convenient to choose a 
"symmetric gauge" given by 

A=\{-xz)xX=^x~^y (58) 

Firstly, we define the kinetic momentum operators as 

p = rjX, so that 11 ~ p + A. Notice the following non 
vanishing commutators 



[P] , Pk] = -ifix^ok j, k e {x, y}. 



(59) 



mality of the Hilbert basis (and its spatial derivatives) we where ejk is the antisymmetric tensor e^y 
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Secondly, we introduce the operators 



1 

2% 
2h^ 



{Py + ip^) 



(Py - iPx) 



(60) 



which satisfy standard commutation relations for Bose 
operators, [a, a^] = 1- The number operator 
Na = a) a satisfies commutation relations [Na,a] — 
—a, [Na, a^] = and we have the identity 



2-q 



(61) 



In analogy with the case of a charged particle in the 
electromagnetic fielcP^, we obtain that the gyrotropic 
translational group is generated by T = 11 + A, 

Tx^Px + xy, Ty = Py - XX (62) 

which satisfies the following commutation relations, 

[Tj,Pk] = 0, [Tj,Tk] = ihxejk, j, k € {x, y} (63) 

Now we introduce another set of Bose operators 

IT 



2hx 



{Ty ~ iT^) 



(64) 



6t = ^— (T,+.T.), 

which satisfy commutation relations 

[b,b^^l, [Ah,b]=-b, [Mb,b^=b\ 

where Mb = b^b is the corresponding number operator. 
Notice that the commutation relations [a, b] — [a, 6^] = 
also hold. 

Coordinates x and y can be expressed in terms of the 
above Bose operators: 



X^^iPy- Ty) 

A 

y {px -T^) 

X 



(65) 



so that 
1 



rjujli (x^ + y^) = hcoM {Na + Mb - ab - a^6^ + l) 



(66) 



Consequently, the Hamiltonian (|57|) becomes 



(1 + ONa + (,Mb - ^{ab + fltfet) +^+\ 



2 
(67) 

It can be diagonalized with the help of Bogoliubov trans- 
formations 



a = ua 



vb\ 



(3 = ub — va' 



with M, V being real numbers. These new operators satisfy 
Bose commutation relations if — = 1. Substituting 



the above equations into Eq. ( 67 1 we obtain 



n' = nuj 



M 



(69) 



(a^^^ + ap) {uv{l + 20 - ?("' + t-')) 
+ {v^{l + 20 - 2^™ + (e + 1/2)) 

To get a Hamiltonian in the oscillator form, the coeffi- 
cient related to (a^/3^-)-a/3) should be zero, which requires 

to(1 -h20 -C(u^ = (70) 

The solution is u = cosh{9), v = sinh(6'), 

tanh(20) = (71) 

Finally, the coefficients of the terms a^a and /3^^ become 

"1 + 2^ 



1 
2 

1 r 

2 
1 
2 

1 r 

2 



cosh(20) 

v/1 + 4e + 1 

l + 2g ^ 
cosh(26l) 

Vl + 4C - 1 



(72) 



and, consequently, the Hamiltonian in the second quan- 
tization formalism becomes 

n' = ni^+(a^a + ]A +hi^_[^p'^p + ]A (73) 



where 



+ 4C ± 1 



(74) 



Noticing that for any n G Z+ we have ^ = £o+ p^^^2 gn, 
the second quantization procedure yields the following 
form of Hamiltonian ( 56 1 



n>0 



n ( ^Ai^n + 



n>0 



(75) 

where oj^ are the eigenfrequencies of the vortex state 
given by 



1 



'(l + 4eo) + 



4A 



Pm^m 



2 9^±1 



(76) 



which coincides with ( 46 ) 
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VI. COMPUTATION OF THE VORTEX MASS 



As it has been discussed in Section II, to calculate the ^i^, — A 



vortex mass density tensor (see Eq. ( 10 1) we need to find 
solution ($(r, t), 0(r, t)) of the Landau-Lifshitz equation 
in the slow dynamics regime, i.e. in the first order on 

\X\. A more convenient set of variables for this problem 

is the pair (<I>,rn), where m = = = cos0 is 

the projection of the magnetic moment onto the z axis. 
Notice that Landau-Lifshitz equation can be recast as the 
set of equations 



~dt 
dm 
~dt 



Ms 5 m 
_J_5£_ 



(77) 



The total energy f ($, m) splits into the sum 



£:($,m) =Eex{^,m) 

Ml 
1 



-£"„„($, m) + £demag{^,m) 

sin^ e(V3>)2] 



1 



= A 



1 — TO" 
,2 



M| 2 
(Vto)2 + (1 - to2)(V$)2 



if„cos^$(l 



K±m~ 



\M-Hd (78) 



with A being the exchange constant, K± being the 
anisotropy constants and being the demagnetizing 
field. Recall that Hd{r) = -V$d(r)> with \/^<i>d{r) = 
—^TTpd and pd = —V • M. Equivalentlj^, 



<S>d{r}= / dVM(f') • V 



d-^r 



,V'-M(r=') 



V 

dV 



dS" 



Al(r') 



3^ Pdir') 



IdV 
dV 



2-y crd(^) 



(79) 



dV 



with ad = M ■ ft being the effective surface "charge" den- 
sity and V being the volume of the system. Consequently, 
the demagnetizing energy can be written as 



1 



dVprf(r)$d(r) 



1 



dV 



d^rad{r)^d{r) 

m 

We are dealing with a two-dimensional micrometric ob- 
ject, so the surface energy term dominates over the vol- 
ume energy term. We can approximate this surface term 
by an effective easy plane anisotropy contribution given 

by 



£deraaa,S^ I dV27rM2(r) 

Jv 



(81) 



This gives for the total energy 
1 



1 



(Vm)^ + (1 - m^){V^Y 



K„ cos^ $(1 - TO^) + {Kj_ + 2TTMt)m^ (82) 



and the equations of motion ( 77 1 become 
Ah d$ 



2Am 



7 dt 



Ms dm 
~~dt 



r(VTO)" 



2A 



A TO 



(1 — m^Y 1 ^ ^^^^ 

- 2Am{V<^f + 2K„ cos^ $ to 
+ 2{K±_ + 2TTM^)m, 

- K„ sin(2$)(l - m^) - AAmVm ■ V$ 
+ 2A(1 -m^) A$ 



(83) 



In the slow dynamics regime {\X\ ^ 1) solutions (<i>,TO) 
can be split into $ — $o + '^'i and to = TOq -I- TOi, 
where $o and toq are the static solutions of the Landau- 
Lifshitz equation (we consider the anisotropy interaction 
to be weak enough so that the static solutions of the 
Hamiltonian E^x + £demag are valid for our problem) and 

where $i and toi are linear on \X\. Static solutions are 



given by Eqs. (24 1. As discussed in Section III, in the 
weak bending regime we can approximate f ~ r so that 

ca, , ^ dm,Q ^ 

v<E>o = — and Vmo = — ; — 



dr 



Linearizing Eqs. (83 1 and taking into account that 

d^ ^ dm ^ 

— ~ —X ■ V$ and that = —X ■ Vto, we obtain the 

. _ dt 

equations of motion 



Ms ^ 60 

UyX ■ ^ ^ - 

7 r 1 



-2 A 



A mi — 



2A ( dmc 



2A 



(1 — TOq)^ \ dr 
2K„ cos^ $0 - 2{K±_ + 2tt M'^s) 



mi 



AA mo dmo ^ „ , . e<<, 

" "i^^^^r ■ ^TOi - AAriy Too — • V'I>i 

(1 — TOg)-^ dr r 

- 2K„mo sin(2$o)<i'i 

Ms ^ ^ dmo 2n 
7 dr 



= 2A{1 - m^o) A $1 



AAriy moVmi 



. , dmo ^ ^ ^ 
AA mQ——er ■ V$i 
dr 



-\- 2K„ sin(2$o)mom,i — 

2A'„cos(2$o)(l - ml)^i 
(84) 

Asymptotic expressions for the 0(|X|) corrections to 
the out-of-plane vortex shape can be determined by sub- 
stituting Eqs. (24) into Eqs. (84). In doing so we obtain 



X-e, 



^ _Ms_ ^ 

27 (if^ + 27rAf|) + cos2 $0 ' 
. C^Ms 3/2.^ . cxp(-r/Ao) 
*^ = ^^° (^-'^^ 



1 



(85) 
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for r ^ Aq ; and 



MsCin^ ^ . 3 
7 18A 



(86) 



for r <^ Aq. Computation of the mass of the vortex core 
can be made via lit, which should be proportional to X 



in this limit: 



7 

Ms 

7 



d^r(V'I>)T 



J dV(V$o)mo J dV(V$o)m: 

y dV(v$i)TOo-^ y" dV(v$i)TO 



Notice that 



7 



d r(V<i>o)TOo = because it corre- 



sponds to the momentum of the static solution. The last 
term of Eq. ( 87 ) can be neglected because it is quadratic 



in \X\. Therefore it remains to calculate the second and 
third terms, which are given by 

[ d2r(V$o)mi = / d2f(Va>o)mi 

1 J 7 Jr<A„ 



271 



7 



/ dV(V$o)™ 

Jr>Ao 

1/4 



A'j_ + 2 ttM J 



ln(i?/Ao) - ^ 1 ^ 



and 



d2f(V$i)mo = : 



27r 



7 

11 2 

504 ^ 49 



/ d^r(V$i)mo 

Jr<Ao 

/ d2f(V$i)TOo 



(89) 



respectively. Notice that ^ 



^Z"^" , exp(-22;) 

da; 

1 X 



°° exp(-2x) 

dx = 0.049 because we are interested in 

1 X 
the Hmit i? > Aq. 

Collecting all terms for the momentum, we get for the 
total mass density 



27r 



M. 



.f2 



ln(i?/Ao) 



0.0014 



(90) 



Notice that we are interested in the limit R ^ Aq, so 
that the term involving ln(i?/Ao) is the dominant one. 
Furthermore, redefining the exchange length by a factor 
close to unity we can always absorb the small numerical 
constant in Eq. (90 1 into the logarithmic term. Magneto- 



crystalline anisotropics, if they are sufficiently large, de- 
stroy the circularly polarized state. Consequently, ma- 
terials like permalloy, used in the studies of the vortex 
state, have negligible magneto-crystalline anisotropy en- 
ergy as compared to the demagnetizing energy. This 
means that the above expression for the vortex mass den- 
sity can be reduced to 



Pm 



1 

47' 



2 ln(i?/Ao) 



(91) 



With account of this formula one obtains the following 
expressions for the parameters ujm and eg that determine 



(87) eigenfrequencies in the equation ( 76 1 



Co 



ln(i?/Ao) 



(92) 



VII. EFFECTS OF THE MAGNETIC FIELD 
AND DISSIPATION 

In this section we study the effects of a magnetic field 
on the excitation modes of the vortex state. Arbitrary 
directed magnetic field can be split into two components, 
one being in the plane of the disk and the other one being 
perpendicular to it. The effects of these two components 
can be investigated separately. 

Consider first the case of a spatially uniform in-plane 
magnetic field, Hin = h^Sx + hyCy. For small displace- 
ments along the disk, the magnetic vortex develops an 
in-plane magnetization density given b}El 

a1{X) = -fi\zx x], fi= {2n/3)Msn^R. (93) 

The Zeeman energy density term is 

LJziX) - -M{X) ■ = -/i [z X i?,;„] • X 

= fihyX - fih^y (94) 
and thus the total in-plane potential energy becomes 



1 



zPmi^m^o 



y 



pK \ 

y 2 — 



2 " 



1 



2 Pm^^m^o 



-Hi, 



(95) 
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Notice that by shifting the origin of the coordinate sys- 
tem we retrieve the original in-plane term of the to- 



1 



tal energy density (39) except for the constant term 
-Hf„, which is field dependent. Conse- 



2 Pm^^jEq 

quently, the application of an in-plane magnetic field does 
not modify the excitation modes given by ( 76 1 . 



Consider now the effect of the magnetic field perpen- 
dicular to the plane of the disk. H± = Hz. Application 
of such a field results in the precession of the magnetic 
moment of the vortex about the direction of the field, 
described by the Landau-Lifshitz equatiorPSl 



dM{t, X) 
dt 



-7 



M{t,X) X Hj_ 



(96) 



where 7 is the electron gyromagnetic ratio. Formally, this 
effect can be accounted for by adding an extra term to 



the gyrovector. Indeed, integration of Eq. (27) (with no 
potential energy) on time gives X = a 
a = -1/pm- 



X X PG 



where 



With account of Eq. (|93|, we have 

X 







zx X 


= -7 


zx X 


a (^z X 


X 


X PG 




= -7 


zx X 



(97) 
(98) 



The vector identity ax b x c— {a ■ cjb— {a ■ b)c leads to 
apG — —jH. Consequently, the precessional effect of the 

perpendicular field can be absorbed into the gyrovector 

^ 7 ^ 
density if one adds to it the term p^ ~ H± ~ 

^ (J,-Hj_ Q, 

PmjH±. This adds the Larmor frequency to ujm- 

LOM (H) = = LOM + = + iH (99) 



Pm 



Pm 



so that the eigenfrequencies (|76| become 
1 



'{l + Ae{H)) + 



4A 



Pm^m{H) 
with e{H) given by 



e{H) 



iOM{H) 
(100) 



lom{H) PMioliiH) (l + 7^M/)2 



(101) 

Introducing dimensionless variables h = 'yH /lom and 
w^(/i) = uj^{H)/ujm we can rewrite Eqs. (100) as 



4eo 



4A 



il ± 1 



(i + h) 



4en 



+ 



sgn(l + h) \ql 
y(l + /i)2+4eo PMU^li 



(102) 



The distance between w+ and oj~ equals Aw = ujm + iH. 

To conclude this Section, we investigate the effects of 
the dissipation on the excitation modes of magnetic vor- 
tices. We consider only the zero field case. Derivation 
of the corresponding expressions when a magnetic field 
is applied is straightforward. The way to introduce dis- 
sipation into our equations is by adding a damping term 

of the fo rm — DX {D being the damping constant) to 
Eq. pTft^^l^si Therefore, the elastic Thiele's equation 

becomes 

PmX - \dlX + Xxpg-dJI + pM^lieoX = (103) 

Repeating the procedure of Sec. IV with the above equa- 
tion in the massive vortex case (pM 7^ 0) we obtain the 
following equation for the frequency modes 

+ {ujm + id)uj - ujlie{q) = (104) 

with d = D/ Pm and e[q) = cq + ^^^'^2 g^- The (complex) 
roots of this equation, uj± = Re(a;±) +ilm(a;-|-), are given 

by 



Re(w±) 
Im(tj±) 



T^cos(0/2) - — , 
d 

T^sin(0/2)-- 



(105) 



with 



[{I + Ae{q))Loli ~ d^Y + Ad^uj 



M 



arg 
arctan 



( [(1 + H<i)Wm - rf'] + * [2dwM] 



2diUM 



{l + Ae{q))ul,-d^ 



(106) 



In the regime of weak dissipation, d << ujm, we have 

and r ~ {1 + 4e(g))wJ,^. As 



arctan 



(1 + 4e{q))LJM 
cos[arctan(a;)/2] ~ 1 
— h o{x^) if |a;| <^ 1, we finally obtain 



X 

~8 



o{x'^) and sin[arctan(x)/2] 



Re(w±) = T 



i(Vl + 4e(g)±l 



UJM 



X 



VI + 4eo 



1 + :r 



3 {d/ujM) 



1 {d/ujM? 
4(l + 4e(g))3/2 

UJM {d/uJMY 

4 (l + 4eo)3/2 

„2 



UJm 



2(1 



Pm^^m 



and 



- 1 



(107) 



(108) 



Im(c^+) ^ 1 + v/l + 4£(g) ^ (1 + y/l + 4e(g))^ 
Im(cj_) 1 - ^l + 4e{q) 4e(q) 

(109) 
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VIII. CONCLUSIONS 

We have studied excitation modes of vortices in cir- 
cularly polarized mesoscopic magnetic disks that corre- 
spond to the string-like gyroscopic waves in the vortex 
core. This problem was studied by classical treatment 
based upon Landau-Lifshitz equation and by quantum 
treatment based upon Hamiltonian approach. The quan- 
tum problem is interesting on its own as it is equiva- 
lent to the problem of quantum oscillations of a charged 
string confined in a parabolic potential and subjected 
to the magnetic field, which in its turn, is a generaliza- 
tion of the problem of the field-induced orbital motion 
of the electron in a potential well. Both treatments ren- 
dered identical results. Our solution generalizes the ex- 
pression for the frequency of the gyroscopic motion of 
the vortex for the case of the finite wave number g, as 
a;_(g) = ujq -|- 7Ms(qAo)^ ln(i?/Ao), where ujg is the 
conventional gyrofrequency, 7 is the gyromagnetic ratio, 
Aig is the saturation magnetization, Aq is the exchange 
length, and R is the radius of the disk. This expression is 
valid in the long- wave limit gAp ^ 1. The wave number 
is quantized, q„ = 2-Kn/ L, where L is the thickness of the 
disk and n is an integer. For a disk of radius R ~ 1/im, 
thickness L ^ lOOnm, exchange length Aq ~ 5nm, and 
saturation magnetization Ms ^ lO'^emu, the n = 1 mode 
is separated from ujg by a few GHz. It could be excited 
by, e.g., a tip of a force microscope or a micro-SQUID 
placed at the center of the disk. Such measurement, while 
challenging, is definitely within experimental reach. 

Throughout this paper we considered of a non-zero 
mass of the vortex. In addition to the gyroscopic mode 
the finite provides a new excitation mode, ^+(9) = 
L0},[ +7Ms(gAo)^ ln(i?/Ao). The gap, wm is higher than 
the gyroscopic frequency ujq. It depends explicitly on 
the vortex mass. The vortex mass density has been com- 



puted by us as a coefficient of proportionality, pM , in the 
kinetic energy of the moving vortex pmv^ /2. It is given 
by pm — 1/(47^) ln(i?/Ao), where R is the radius of the 
disk, Aq is the exchange length, and 7 is the gyromag- 
netic ratio. For a 25nm thick, micron size permalloy disk 
this gives the vortex mass in the ball park of lO^^^kg, 
which is close to the experimental value estimated for 
a comparable size permalloy ring^^. Our result for the 
mass gives ujm — Sir jMs / ln{R/ Aq) . This is in the ball- 
park of, or below, the uniform ferromagnetic resonance 
of the disk. It would be interesting to investigate this 
frequency range experimentally alongside with the low- 
frequency gyroscopic mode. One can also test in experi- 
ment the explicit field dependence of the vortex modes, 
computed in this paper. So far we have done it for the low 
field that only slightly disturbs the vortex state formed 
in a zero field. However, the statement concerning the 
existence of the additional mode due to the finite vortex 
mass should apply to higher fields as well. This case, 
however, defies analytical study and must employ full- 
scale numerical micromagnetic calculations. When the 
field is sufficient to fully polarize the disk in the perpen- 
dicular direction, we expect the high frequency mode to 
evolve into the uniform ferromagnetic resonance. 
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